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Abstract 

The observed CMBR dipole is generally interpreted as a Doppler 
effect arising from the motion of the Earth relative to the CMBR 
frame. An alternative interpretation, proposed in the last years, is that 
the dipole results from ultra-large scale isocurvature perturbations. 
We examine this idea in the context of open cosmologies and show 
that the isocurvature interpretation is not valid in an open universe, 
unless it is extremely close to a flat universe, |flo — 1| < 1CP 4 . 


1 Introduction 

In the standard Big Bang scenario, the oldest relic from the early universe 
is the Cosmic Microwave Background Radiation (CMBR). This radiation is 
remarkably close to isotropy. When it is decomposed into multipoles, 
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one finds that the dipole contribution is of the order (A T/T)i=% ~ 10 -3 , 
while the contribution from higher multipoles is only (AT/T)oi ~ ICR 5 (see 
e.g.|Ij]). Usually, the dipole term is interpreted as a Doppler effect, i.e. as 
the consequence of our local motion with respect to the “CMBR rest frame”, 
whereas the other multipoles are accounted for by primordial cosmological 
perturbations. However, one cannot reject a priori the possibility that a sig¬ 
nificant part of the dipole originates as well from cosmological perturbations. 
In fact, this alternative idea was stimulated by the results of Lauer and Post¬ 
man [j|| in 1991, who found a dipole in distant Abell clusters inconsistent 
with the CMBR dipole. Other observations based on nearby galaxies, IRAS 
galaxies or distant supernovae tend on the other hand to favour the orthodox 
interpretation, but are still inconclusive. 

Let us be more precise on the distinction between a cosmological dipole 
and a local dipole. If our universe was perfectly homogeneous and isotropic 
then there would be no cosmological dipole; a comoving observer would not 
detect any dipole whereas a non comoving observer would measure a dipole, 
in this case a purely local dipole, due to his motion. In a perturbed model, 
it is more delicate to distinguish between a local and a cosmological dipole. 
The perturbation of any quantity can be decomposed into contributions from 
different scales. The cosmological dipole is then due to very large scale 
perturbations, typically scales of the order of the Hubble radius, and the 
local dipole is due to the small scale component (the dominant one) of the 
observer peculiar velocity. 

The purpose of this paper is to explore the origin of the dipole from a 
theoretical point of view and to ask whether it is possible, in the context 
of Friedman-Lemaitre-Robertson-Walker (FLRW) cosmologies with gaussian 
random fields of linear perturbations, to obtain a cosmological dipole and 
higher multipoles compatible with the observations. To discriminate be¬ 
tween a local Doppler dipole and a cosmological dipole is more, in our view, 
than a mere academic exercise because the Doppler assumption enters in the 
quadrupolc analysis of the observations. The reason is that our motion would 
not only induce a dipole but also a quadrupole, and other multipoles. 

It was shown by Paczynski and Piran [[| that a strong dipole can be found 
within a particular Tolman-Bondi model with an inhomogeneous radiative 
component. In fact, as was shown recently by Langlois and Piran JiJ, this 
result can be obtained, more generally, in a perturbed FLRW flat cosmology, 
from ultra large scale (larger than the current Hubble radius) isocurvature 


2 


perturbations. Here, we reinvestigate this possibility in the context of open 
cosmologies. 

The influence of stochastic cosmological perturbations with wavelengths 
larger than the current Hubble radius on the CMBR temperature fluctuations 
was examined initially by Grishchuk and Zel’dovich ||. Their analysis was 
restricted to perturbations in a flat background for which the Fourier expan¬ 
sion is available. When one considers an open universe one needs the com¬ 
plicated formalism of mode functions in an hyperbolic space. Recently, Lyth 
and Woszczyna || stressed the importance of super curvature modes, which 
had been ignored previously, when dealing with random fields. With this 
new element, the influence of ultra large scale perturbations on the CMBR 
for open universes was then investigated in references [fj and [T|. The present 
work extends these papers in two respects. First, we focus our attention on 
the dipole, which was not treated previously, and, second, we allow not only 
for adiabatic perturbations but also for isocurvature perturbations. 

The paper is organized as follows. In section 2, we recall the main steps 
of the derivation of the Sachs-Wolfe effect, responsible of the large angle 
anisotropies in the CMBR. Our approach follows the work of Panek [§], 
although our notation is different. The final result is valid for flat or curved 
FLRW background geometries and is expressed in terms of gauge-invariant 
variables. Section 3 gives the CMBR multipoles in a flat geometry. Although 
these results were already given in [|], the actual presentation is original 
and provides a bridge to the formalism needed for an open cosmology. In 
section 4, the case of an open cosmology is considered. Our notation follows 
essentially Refs || and Finally, section 5 contains the conclusions of this 
work. 


2 The Sachs-Wolfe effect 

Our background geometry is described by a FLRW metric, 

ds 2 = a 2 (if) dij 2 + n / i jdx l dx j ') , ( 2 ) 

where dl 2 = 7 ijdx l dx J is the metric of a flat or curved maximally symmetric 
space (a flat space corresponds to K = 0 and an open space to K = — 1, 
where K is the normalized space curvature). The CMBR is composed of 
cosmological photons, which became free after the recombination of protons 
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and electrons and the subsequent decoupling of matter and radiation at a 
redshift zi s . For numerical applications we shall take zi s = 1000 (zi s depends 
in fact on the spatial curvature but this dependence is very weak [j9J and 
can be safely ignored here). All the photons reaching us now were emitted 
on a physical surface at the epoch of decoupling, called the last scattering 
surface, and defined physically by n e = const , where n e is the density of free 
electrons. The universe will be supposed to be matter dominated from the 
last scattering time to now, which is an excellent approximation except for 
low values of h2 0 . 

Now, we consider a perturbed FLRW universe, endowed with the metric 
(in the longitudinal gauge) 


ds 2 = a 2 (r/) —(1 + 2 ^)drj 2 + (1 — 2T)7 ijdx l dx j 


(3) 


Only perturbations of the scalar type will be considered here (see e.g. E3) 
and, for simplicity, the anisotropic stress of the matter will be supposed to 
vanish so that the Einstein equations yield 


<3> = T. 


(4) 


Because of the presence of geometric perturbations, photons will be redshifted 
with slight differences depending on their position of emission on the last 
scattering surface. The resulting fluctuations in the temperature measured 
by an observer were first calculated by Sachs and Wolfe |nj. The origin of 


this effect is purely geometric and is dominant only for large angular scales, 
which correspond essentially to scales larger than the Hubble radius at the 
epoch of last scattering. On smaller scales, causal processes must be taken 
into account. In this work, we are interested only by the first multipoles of the 
CMBR, which are dominated by the last scattering superhorizon modes, so 
that only the Sachs-Wolfe effect is relevant. In the following, the Sachs-Wolfe 
effect is rederived shortly following the more modern calculation of Panek 
[Sj]. However our notations are different and the final result is expressed in 
another form. 

Denoting the vector tangent to a null geodesic, the evolution of the 
photon is governed by the geodesic equation: 


k v d v k? + r % T k a k T = 0 , 


(5) 
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where the T^ r are the Christoffel symbols. We write the perturbed tangent 
vector k M = {z/(l — M),—ue l + P 1 }, where the terms in v represent the 
unperturbed solution; e l is a unit spatial vector, i.e. such that = 1 ; 

M is the perturbation of the frequency, P l the perturbation of the spatial 
direction. Here, we are interested only in the frequency perturbation and 
thus keep only the time component of the geodesic equation. After using the 
identity g fiu k^k l/ = 0 at first order, which reads 

7 ij e i P j = (M — <L — ) v, ( 6 ) 


in order to eliminate the P l in the time component of the geodesic equation 
at first order, one finds eventually 


dM 

~dX 


(dr, - e%)M = $' - 2e i d& - 'L'. 


(7) 


The prime stands for a derivative with respect to the conformal time r/. Note 
that the parameter A is related to the affine parameter r of the unperturbed 
geodesic by the relation dr/dX — \/v. The prime stands for derivative with 
respect to the conformal time rj. Equation (@) corresponds to Eq. (29) of [§]. 

Typically one can model the content of the universe from the epoch of 
the last scattering to now by two fluids, a pressureless baryon component 
and a radiation component, described locally by a black-body spectrum. Let 
us define the temperature of the photons as seen by observers moving with 
the baryon component, so that the temperature ratio between emission and 
reception is given by 

T r = (k^vJpR 

T e (k^) E [ } 

In the perturbed geometry ([]), one thus obtains 


Tr 

T e 
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where the velocity potential V is defined by 5u\ = adiV. As a consequence, 
the temperature fluctuations on the sky measured by an observer today (de¬ 
noted by the subscript 0 ) are given by 


ST 
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The left hand side term is a function of the celestial coordinates correspond¬ 
ing to the direction of observation for the observer. The right hand side is 
expressed as a function of the emission point, of coordinates (r/ e ,x*), defined 
as the intersection of the last scattering surface with the null geodesic go¬ 
ing through the observer with the given direction , as well as the observer 
position (■ rj 0 ,x l 0 ). In general the physical last scattering surface is distinct 
from the constant time hypersurface rj = r)i s , but one can replace the first 
two terms on the right hand side of ( |10| ) by (ST/T)(rji s , x\ s — x l e ) (i.e. as a 
function of the intersection point between the light geodesic and the hyper¬ 
surface r] = T]i s ). This follows from the local conservation law aT = const 
for free radiation. Then, expressing the term [$ — M]° e in ( jldp as an integral 
over the null geodesic and using equation (j7j), one finds 



= ~ A ^\ ls + (hV) u + [-$ + e%V 
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where h = a'/a. Using the Stefan law p 7 oc T 4 , we have replaced the tem¬ 
perature fluctuations by fluctuations of the radiation energy density and also 
introduced the comoving energy density perturbation 5, which can be de¬ 
rived from the energy density perturbation in the longitudinal gauge by the 
expression S = (Sp/ p)l + (p'/p)V. To obtain (|lT|) , we have used implicitly 
the equality of the velocities of radiation and baryonic matter at the time 
of last scattering, which results from the preexisting tight coupling between 
the two fluids. The term <5 7 /4 in ( |TT| ) represents the intrinsic temperature 
fluctuations on the last scattering surface. The rest of the expression is the 
Sachs-Wolfe effect. In (H]) all the quantities at the last scattering epoch are 
now evaluated on the hypersurface rji s which does not necessarily coincides 
with the physical last scattering surface (but all the corrections would be 
second order perturbations). 

It is instructive to decompose the Sachs-Wolfe expression into several 
components and to rewrite (0)> after dropping the term —$0 which con¬ 
tributes only to the monopole, in the form 


ST 

~T 


R 




( 12 ) 
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with 



int 





= (i> + hv) ls , — 


pSW 



Dop 


and 



isw 


[° dX ($' + T') 

Jls 


= e i d i (V 0 -V ls ), 
(13) 


(14) 


The term ($ + hV)i s will be called the proper (or ordinary) Sachs-Wolfe 
effect (pSW), because only this term was computed in the pioneering paper 
of Sachs and Wolfe 0- The term e l d t (Vo — Vi s ) is the difference between 


the observer velocity and the emission velocity along the line of sight and 
corresponds to a Doppler effect (Dop). Note that the local dipole is due 
only to a term of the form e l v®, where v® is the contribution of small scales 
to our peculiar velocity; this contribution, which comes from the nonlinear 
evolution of the perturbations, will be ignored in the rest of the paper where 
we consider only the effect of very large scales. Finally the integral term is 
usually called the integrated Sachs-Wolfe effect (ISW). 

The linearized Einstein equations and fluid equations (see e.g. Hi 


or 


g) can then be used to express all the terms in (|TT| ) in terms of only 
one quantity, the most convenient being the gravitational potential $. For 
example, the velocity potential V can be related to $ via the Euler equation 
and the conservation equation: 


V = 


2 

3 (h 2 + K) 


($' + M>). 


(15) 


Via the Einstein equations, one gets a relativistic generalization of Poisson 
equation relating the gravitational potential to the total comoving energy 
density St, which reads: 


(A+ 3 K)$ = ^{h 2 + K)8 t . (16) 

At this stage one must distinguish two kinds of perturbations: adiabatic and 
isocurvature perturbations. Any general perturbation can be seen as a linear 
combination of these two kinds. An isocurvature perturbation corresponds 
to a primordial (in the radiation era) perturbation of the matter which does 
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not affect the geometry: it implies a perturbation in the relative composition 
of the cosmological fluid with no perturbation in the total energy density. 
On the contrary an adiabatic perturbation is a primordial perturbation in 
the total energy density without modification of the relative quantities of the 
fluids. For an adiabatic perturbation, the relation between the matter and 
radiation density perturbations is by definition <5 7 = (4/3 )S m , which implies, 
during matter domination, 

S-r - ( 17 ) 

As will be seen later, the consequence of this relation together with the Pois¬ 
son equation (]T6|) is that, both for the flat and open cases, the intrinsic fluc¬ 
tuations resulting from adiabatic initial perturbations on scales larger than 
the Hubble radius are always negligible. For an isocurvature perturbation, 
during matter domination, 

- ~S, ( 18 ) 

where S — S m — (3/4)S 7 is the entropy perturbation. Both for open and 
flat universes, the appropriate combination of the equations of motion (i.e. 
conservation and Euler equations) of the two fluids, radiation and baryons, 
tells us that S is constant on scales bigger than the Hubble radius, which is 
the case for perturbations we are interested in. Isocurvature perturbations 
also produce a gravitational potential. Using the conservation and Euler 
equations for the two fluids, one obtains that the gravitational potential is 
given, at the beginning of the matter era, by 

$ = (19) 

5 

for the modes larger than the Hubble radius. 

To summarize, the multipole coefficients of (P can be decomposed into 
the sum of contributions 


( 20 ) 


The rest of the paper will consist in evaluating the various terms of this 
expression for open/flat universes with adiabatic/isocurvature perturbations. 
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3 Flat universe 


In this section, we consider the flat FLRW models. To study the CMBR 
fluctuations, it is more convenient to work in spherical coordinates with the 
observer at the center, naturally adapted to observations on the celestial 
sphere, rather than the more usual cartesian coordinates: 

dl 2 = dr 2 + r 2 (d6 2 + sin 2 Odcp 2 ') . (21) 

Any perturbation /(r, 9, <fl) can then be decomposed into 

f(r, 6 , (p , rj) = f dk E f klm(jl)Q klrm (22) 

Im 

where the mode functions Qkim are the eigenfunctions of the spatial Laplacian 
with the eigenvalues — k 2 and are defined by 

Qkim = ]j^kji(kr)Y lm (9,<f)), (23) 

where the j) are the spherical Bessel functions and Yi m are the spherical 
harmonics. The definition of the Qkim is such that they are normalized i.e. 

J r 2 dr sin 6 d6 d<j)Q* klm Qk>v m > = 5(k - k')Su'S mm >. (24) 

This decomposition is the equivalent of the Fourier decomposition but in a 
system of spherical coordinates. This presentation is a good introduction to 
the open case where the spherical coordinates are the most natural spatial 
cooordinates. 

All the quantities introduced in the previous section can now be decom¬ 
posed on the mode functions Qkim■ As announced before, we shall focus 
essentially on the gravitational potential <F. The linearized Einstein equa¬ 
tions yield a second order differential equation for the evolution of each mode 
Qkim, which shows that the modes larger than the Hubble radius, i.e. k < h, 
are constant in time in a flat universe during the matter era. This result has 
several consequences. The first is that the ISW contribution vanishes. The 
second is that the relation (pA|) simplifies to 

hV = - 2 -<$> 1 (K = 0). (25) 
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Now, using the orthonormality property ( |24|) with equations (|TT| ) and (|25|), 
the Sachs-Wolfe harmonic coefficients can be expressed in the form 

<c = /! r dk ^yvr(k), (26) 

with the window function Wj sw (k) = Wf SU ( k ) + W[ >op (k) 1 the proper SW 
and Doppler contributions being respectively given by 

Wf SH (k) = ji(kr ls ), W° op (k) = 2hJ s 1 -^-j l (kri s ) - 2Hq l ^ji(kr 0 ). (27) 

ri s is the comoving distance from the hypersurface t] = t]i s . It is related to 
the redshift zi s (= ( a 0 /ai s ) — 1) by the relation 

r is = 2 hg 1 (l - (1 + z ls y 1/2 ) = 2 (Tig 1 - h^ 1 ) . (28) 

Using the relations 

4 ~ji(kr ) = kji-i(kr) - l -^-ji(kr), (21 + l)ji(x) = x(j i+1 + jz- i(x)), (29) 

dr r 

one can rewrite the dipole term 

W? w (k) = -khs 1 (jo(kr ls ) - 1) + ^(1 - 4htt)j 2 (kn.). (30) 

and the terms l > 1 

/ 9/ \ 2k 

Wi W {k) = ( 1 + ) Ji( kr is) ~ -j^ji+i(kris). (31) 

Since jo(x) = sinx/x and ji(x) ~ x l for x « 1, the dipole term Wf w (k) be¬ 
haves like (kri s ) 3 for small k. On the other hand, for l > 1, Wf" ( k ) goes like 
(kris) 1 for long wavelengths. This means that the Sachs-Wolfe dipole on ultra 
large scales is suppressed with respect to the quadrupole. Let us now consider 
the intrinsic contribution. In the case of adiabatic perturbations, equation 
(1171) and the Poisson equation (pTJ) (with K = 0) shows that the intrinsic 
contribution is always negligible with respect to the Sachs-Wolfe contribution 
for modes larger than the Hubble radius at the time of last scattering. For 
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adiabatic perturbations, the temperature fluctuations reduce to the Sachs- 
Wolfe fluctuations. The (total) cosmological dipole is thus suppressed with 
respect to the quadrupole (this remarkable property was stressed in |T^]). If 
this case describes the reality we live in, the conclusion is that the observed 
dipole is necessarily of local origin. 

For isocurvature perturbations, the intrinsic contribution cannot be ne¬ 
glected and the total window function Wi(k) is the sum of the SW window 
function and of the intrinsic window function which, using (O) and 


as 


well as ([191) , is found to be 


Wi so (k) = 5 j t (kr ls ). 


(32) 


Hence Wi behaves like kri s for small k and can be made as large as one 
wishes with respect to the quadrupole window function by considering large 
enough scales. 

The above analysis has been dealing with individual modes. To get the 
total multipoles, one must sum on all the modes. Assuming that $ is a 
gaussian random field described by the power spectrum 

i^kim^k'Vm') = 27r 2 fc 3 Vq>(k)S(k — k')5a>b mm i (33) 


the expectation value of the various multipoles reads 

4-7T r die 

(k™! 2 ) = -j ] (34) 

Once the power spectrum is given, the difference between the multipoles 
comes only from the window functions. The most common power spectrum 
is the flat spectrum corresponding to a constant V<i>(k). On Fig. 1, we 
have plotted the square of the dipole and quadrupole window functions in 
the adiabatic case. Since the dipole window function is always subdomi¬ 
nant with respect to the quadrupole one or at best of the same order of 
magnitude for k ~ h 0 , the conclusion is that the observed dipole can only 
be explained by a local effect in that case. On the contrary, for isocurva¬ 
ture perturbations, Fig.2 shows that the dipole window function dominates 
on ultra large scales and the observed dipole can be of cosmological origin. 
However the corresponding isocurvature spectrum cannot be constant on all 
scales because the main contribution to the dipole and quadrupole would 
come from scales k ~ h 0 and the resulting dipole and quadrupole would be 
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Figure 1: squared amplitude of the window functions in the adiabatic case 
for the dipole (solid curve) and the quadrupole (dashed curve). On ultra 
large scales the dipole window function is suppressed with respect to the 
quadrupole one. 



k/hO 


Figure 2: squared amplitude of the window functions in the isocurvature 
case for the dipole (solid curve) and the quadrupole (dashed curve). On 
ultra large scales the quadrupole window function is suppressed with respect 
to the dipole one. 
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comparable. To obtain cosmological dipole and quadrupole compatible with 
the observations, two conditions must therefore be satisfied: the existence of 
ultra large scale isocurvature perturbations; and their spectrum must be sup¬ 
pressed for scales A < lOOff^ -1 . Although this cut-off condition can appear 
on first thought difficult to fulfill in a sensible model, it turns out that two 
relatively simple ideas can produce such an effect. One idea is the existence of 
pre-inflationnary isocurvature perturbations, which would have been pushed 
away far beyond the horizon during the inflation era jnj; the second idea 


is to consider the simplest two scalar field inflation model, which naturally 
produces an isocurvature spectrum with an abrupt cut-off 
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4 Open universe 

We now write the background spatial metric in the form 

dl 2 = dr 2 + sinh 2 r(d0 2 + sin 2 ddcj) 2 ). (35) 


With this choice of coordinates the Friedmann equation reads 


H 2 


8nG 



(36) 


where H = ah is the Hubble parameter. As is clear from the above formula, 
a represents the curvature scale: on scales smaller than a, the curvature of 
space is not “felt” whereas scales larger than a are affected by the curvature. 
It is usual to parametrize the space curvature by the ratio of the energy 
density with respect to the critical energy density today, Oo = 3Hq/8ttG. 
Oo = 1 for a flat space and fio < 1 for an open space. It then follows, in the 
matter era, that the comoving Hubble parameter can be expressed in terms 
of the redshift 0 and 


h 


1 -f- 

1 — f2 0 


(37) 


In particular h 0 = 1/y/l — Ho represents the ratio between the current curva¬ 
ture scale do and the Hubble scale Hq 1 . For an open universe the curvature 
scale is thus always larger than the Hubble scale and all the more larger as 
fl 0 is closer to one. 
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In an open universe, the Fourier treatment does not apply. One must 
then use the formalism of mode functions for an hyperbolic space. Following 
[6i], any gaussian random held can be expanded in the form 

OO l 

dk E E fkl m {j])Qklmi 

Z=0 m=~l 

where the Qkim are eigenfunctions of the spatial Laplacian with eigenvalues 
— k 2 . They can be chosen of the form 

Qkim — 0 ), ( 39 ) 

where Yi m (6,<f>) are the usual spherical harmonics on the two-sphere. The 
radial mode functions II^z can be classified into two categories depending on 
the sign of q 2 = k 2 — 1. The radial functions for the subcurvature modes, 
corresponding to the values q 2 > 0, are given by 



roo 

Jo 


n H = N M u kh 


(40) 


with 



■ i -I —V 2 

IE « 2 + 9 2 ) 

:n =0 


n ki = q 2 (sinhr) / 


—1 d 
sinli r dr 


i +1 


cos(gr). 


(41) 

The subcurvature modes provide a complete orthonormal basis for square 
integrable functions and for this reason were the only modes taken into ac¬ 
count by cosmologists. However, as stressed recently by Lyth and Woszczyna 
[|], they are not enough to describe properly gaussian random fields and one 
must also include the super curvature modes, corresponding to the values 
-1 < q 2 < 0. Their radial functions are obtained from the subcurvature 
modes by analytic continuation: 


klfcz = N k iTi kh 


(42) 


with 


N k0 — 



N kl = 



ri(« 2 +<t) 


(OO) 


(43) 
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and 


II kl = \q\ (sinhr) 


—1 d 


Z+l 


cosh(|g|r) 


sinh r dr 

The normalization of the mode functions is such that 

n kl (r)Tl kr (r) sinh 2 rdr = 5(q - q')5 U '. 


(44) 


(45) 


-'0 

Moreover for a homogeneous gaussian random field, one defines the power 
spectrum by 


2t r 2 

( fklmfl'Vm /) = ~ k ') 5 H' 5 mm'■ ( 46 ) 

Note that in the flat space limit, corresponding to the limits k —> oo and 
r —> 0 with kr fixed, 

n w (r) -> ^~kji(kr ), ( 47 ) 

and one recovers equation (^) for the power spectrum. 

In contrast with the flat case, the gravitational potential modes <&kim on 
scales larger than the Hubble radius evolve with time. Their evolution is 
given by the linearized Einstein equations (see e.g. [ |T^] ): 

$k(ri) = F(v)®k, (48) 


with 


F(rj) = 5 


sinh 2 rj — 3rj sinh q + 4 cosh rj — 4 
(cosh;/ — l) 3 


(49) 


(with the normalization a = cosh?/ — 1 for the scale factor). F is normalized 
so that F —» 1 when 77 —> 0. is constant in time and is given by the 
initial spectrum at the beginning of the matter era. At that initial stage, the 
curvature can be ignored. In fact, even for small values of Oo, F(rj) remains 
almost constant during most of the matter era and drops suddenly in the 
recent past. Similarly, the velocity potential obeys to a law different from 
the flat case and given by 


hV{rj) 


-fcwi, 


(50) 
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with 


G(v) 


15 (2 ?] + rj cosh r] — 3 sinh rj) sinh rj 
16 sinh 6 ( 77 / 2 ) 


(51) 


Like F, G has been normalized such that G —> 1 when 77 —^ 0 and G is 
significantly different from one only for small redshifts. 

The decomposition of the Sachs-Wolfe terms in (|TTD on the basis Qkim, 
like in (|58D, yields for the harmonic coefficients the expression 


a, m = M- [*° dk\q\li tlm \wtf w + WS°> +Wl? W 

7T JO 


with 


,pSW 

4/ 

,Dop _ r, 7-1 


Wir = (3F la -2G la )Il kl (r la )/\q\, 

2 . 


and 


Wm W = 


W kl ° P = 2h^G ls d r U k i(r ls )/\q\ - Zh^Gokd^ 


-6F ls Il kl (r ls ) + 6 / dr F(r] ls + r ls - r)9 r II w (r) 

JO 


(62) 

(53) 

(54) 

/kl- (65) 


Except for the flat case limit, we shall take Fi s = Gi s = 1 in the above 
expressions, which is an excellent approximation because Zi s » 1. A hat 
means a division by ^2/tt. h\ s is the comoving Hubble parameter at the time 
of last scattering. r\ s is the radius coordinate corresponding to the emission 
of the photons, i.e. the coordinate of the last scattering surface. This radius 
can be related to the redshift £i s via the relation 


sinhrfA) = —(1 + z ) 
ho 


-ln-2 




^ 0 ^ T (^0 — 2) f(l + — 1 


(56) 


We must also consider the intrinsic contribution af^ to the harmonic 
coefficients. For adiabatic fluctuations, like in the flat case, the intrinsic 
contribution is negligible as can be seen from ( pT| ) : the reason is not this time 
the smallness of k 2 but the fact that h 2 >> 1 (see (|37p). For isocurvature 
perturbations, the intrinsic window function is given by 


Wir = 5F ls U kl (r ls )/\q\, (57) 
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Finally, once a power spectrum is given, the expectation values for the 
harmonic coefficients can then be written in the form 

= (58) 
where Wi(k) is the sum of all the relevant window functions. 


4.1 Ultra large scale limit 

In the flat case, it was possible to obtain a dipole much larger than the 
quadrupole with isocurvature perturbations on scales much larger than the 
Hubble radius. We wish now to consider such scales in an open universe. In 
the ultra large scale limit, corresponding to k —> 0, the general expressions 
can be simplified because n fc /(r) ~ kfiir). Following [0], we define 


Nl = k™ kNkl = \n n > 2 - X ) 


- 1/2 


(59) 


for l > 1 and N\ = 2/%. Similarly, we define 

III = lim tl k i/k 2 . 

k —>0 

The expressions we are mainly interested in are 

1 / r 

hli(r) = - cothr- ^ 

H ’ 2 V sinh 2 r 


and 


n 2 (r) = - 


1 + 3 


1 — r coth r 
sinh 2 r 


The window functions ( |53H55[ ) can then be expressed in the limit k 
the form 

Wki ~ kCiiQ, o, zi s ), 

with 

cf w = JV,n, (r„), = 2N,l^ 1 ^-n,(r h ) - hi, l'G 0 S u , 


(60) 

(61) 

(62) 
0 in 

(63) 

(64) 
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Sachs-Wolfe dipole, limit k->0 



omega 


Figure 3: amplitude of the Sachs Wolfe dipole window functions in the limit 
k —> 0 as a function of Q 0 - The total Sachs-Wolfe window function (solid 
curve) is the sum of the pure Sachs-Wolfe contribution (dotted curve), of the 
Doppler contribution (dashed curve), which almost cancel each other, and 
finally of the ISW contribution (dashed-dotted curve). 
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Sachs—Wolfe quadrupole, limit k—>0 



0.2 0.4 0.6 0.8 1.0 

omega 


Figure 4: amplitude of the Sachs Wolfe quadrupole window functions in the 
limit k —> 0 as a function of f2 0 - In the total Sachs-Wolfe window function 
(solid curve), The Doppler term (dashed curve) is very small with respect 
to the two other terms, the pure SW term (dotted curve) and the ISW 
term (dashed-dotted curve). One notices that the SW dipole vanishes near 
fl 0 - 0.3. 
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and 

ms d ~ 

£ ISU = N } ‘-611/(ris) + 6 / dr F(rft s + r ls — r)—II/(r) . (65) 

Jo dr 

These three quantities, as well as their sum corresponding to the total Sachs- 
Wolfc effect are plotted as functions of Qq for the dipole (Fig. 3) and for the 
quadrnpole (Fig. 4). The first striking difference with the flat case is that 
the Sachs-Wolfe dipole is not suppressed with respect to the quadrupole. 
One can see more precisely on Fig. 3 that this is due principally to the ISW 
contribution, even if the pure SW and Doppler contributions do not cancel 
each other exactly, especially for low Q 0 . As for the quadrupole, one can see 
that the Doppler contribution is extremely small in comparison with the other 
terms. Another remark is that the quadrupole vanishes around O 0 — 0.3 for 
which the pure SW effect and the ISW effect happen to compensate each 
other. 

In Fig. 5 and Fig. 6, we have plotted, still in the ultra large scale limit, 
the total window functions for the multipoles l = 1,2,3 with, respectively, 
adiabatic and isocurvature initial conditions. We have added the multipole 
l = 3 to show that even if the dipole can be made much bigger than the 
quadrupole in the zone Do — 0.3, this cannot explain the observed CMBR 
since the next multipole is roughly one tenth of the dipole. The conclusion is 
therefore that, in an open universe, the cosmological dipole is typically of the 
same order of magnitude as the higher multipoles. Within our assumptions, 
the observed dipole can be explained cosmologically only in a flat space. 
In the next section, we examine in detail this flat space limit in order to 
put a limit on the flatness required to explain the dipole by cosmological 
perturbations. 

4.2 Flat space limit 

We saw in the previous subsection that, for supercurvature modes, the win¬ 
dow functions go to zero when D —> 1. We need to examine now in more 
details what is the relative importance of the dipole with respect to the 
quadrupole in this limit. To do so, it is convenient to introduce the para- 
mater 
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ultra large scale limit, 1=1,2,3 (adiabatic) 



Figure 5: squared amplitude of the total window functions (intrinsic and 
Sachs-Wolfe) for ultra large scale adiabatic perturbations (limit k —> 0), for 
the three multipoles l = 1, 2, 3 (from top to bottom). 
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ultra large scale limit, 1=1,2,3 (isocurvature) 



Figure 6: squared amplitude of the total window functions (intrinsic and 
Sachs-Wolfe) for ultra large scale isocurvature perturbations (limit k —> 0), 
for the three multipoles / = 1, 2, 3 (from top to bottom). 
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One can immediately notice that the curvature radius reads a 0 = Hq 1 /e so 
that the curvature radius goes to infinity with respect to the Hubble radius, 
as one can expect in a flat space limit. Expanding (|55|) in e, one finds 


r(fi 0 ,z) = 2 1 


-U i- 


\/l + z / 
1 


e + 

\ 3' 


vTT 


1 + z 
e 3 + 0(e 5 ) 


z + 3 — 3\/l + z + 


2^TT 


(67) 


The expansion of the terms appearing in the total Sachs-Wolfe effect is 




= ~k i 

3 l 


a/T+~Z, 


■-(3 + k 2 ) 1 


2, 


VT+z 


, 11 
(, + T 


22 + 35 ^ 


14a/t+k 


)/(i+*) 


e 3 + 0(e 5 ), 


( 68 ) 


w d° p _ 


Y kl 



4 + 3fcp (1 - 


a/1 + 


'1 1 f 7z-2 

7 + Trn \ 14 (! + z) 


' + 0(e 



(69) 


W 1 ™ = 


16 k 


x/rr 


i + 


\/l + z 


+ 0(e 5 ). 


(70) 


Note that the strict limit Ho = 1 is not obtained simply by setting e = 0, 
which would give zero for all the terms above. Indeed it makes more sense 
to label the perturbation scales by the physical quantity k/h 0 instead of the 
unphysical wavenumber k. Since h 0 = 1/e, the proper flat limit is given 
by e —» 0 with ke fixed. However, this limit applies only to subcurvature 
modes whereas supercurvature modes correspond to larger and larger scales 
as f2 0 —> 1 since they are bounded by k < 1. The most striking feature of 
(|68|-[70j) is that the dominant terms of the pure SW contribution and of the 
Doppler term just cancel each other, ft is thus clear that the suppression 
of the dipole is really a consequence of the flatness. We have also written 
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explicitly the terms at the order e 3 because they represent the leading order 
in the total SW window function. For the quadrupole, one obtains 


= ±k(k* + sr (1 - ^=L =) 2 £ 2 + o(£ 4 ), 


(71) 




—+ 3) 1/2 ( yj=j=j - 

W’J W = 0(e 4 ). 


1+z, 


+ 0(£ 4 ), 


(72) 

(73) 


For isocurvature, one must add the intrinsic contribution which is roughly 
five times VV|?f W . By comparing (|68|- ff0|) with ( |7Tl - f73l) , one can see that for 
k « 1 all the window functions, both for the dipole and the quadrupole, 
are proportional to k. As a consequence one cannot obtain a significant 
dipole by considering sufficiently large scales, like in the flat case. The only 
possible difference between the dipole and the quadrupole can thus come 
only from the parameter e. For adiabatic perturbations, the cancellation 
of the e terms implies that the dipole is e times the quadrupole. On the 
contrary, for isocurvature perturbations, the total dipole is proportional to e 
and can be much bigger than the quadrupole which goes like e 2 . To obtain 
a dipole 10 2 times larger than the quadrupole, one thus needs e < 1CT 2 , i.e. 

|ft 0 -1| < io~ 4 . 


5 Conclusions 

The main conclusion of this paper is that a cosmological origin for the dipole 
(or a significant part of it) must be rejected for an open cosmology, which 
means that the observed dipole would be in that case essentially a Doppler 
effect resulting from our local motion dominated by small scale density per¬ 
turbations. This conclusion assumes that the universe can be described as 
a FLRW model with gaussian fields of linear perturbations. More exotic 
possibilities have not been considered here. 

Another conclusion of this work is that the suppression in a flat universe 
of the dipole due to ultra large scale adiabatic perturbations, with respect to 
the quadrupole, is no longer true in an open universe. 
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It turns out that the only possibility to get a dipole two orders of mag¬ 
nitude bigger than the other multipoles is a model with isocurvature per¬ 
turbations on scales of the order of 10 2 Hq 1 in an almost flat geometry: the 
constraint is (Q 0 — 1| < 10 -4 . This constraint can also be interpreted as 
demanding that the curvature scale be one hundred timed larger than the 
Hubble radius. When this constraint is satisfied but with fl 0 11011 strictly 
equal to one (this would be the case for example of an open universe after 
a phase of inflation), two kinds of perturbations can produce the observed 
dipole : the subcurvature modes on scales larger than one hundred times the 
Hubble radius, like in the flat case; but also the super curvature modes, which 
do not exist in the flat case. 

To conclude, a reliable measurement of h2 0 ( see 101 f° r a discussion on 


measuring h2 0 ) different from one would strongly confirm the Doppler assump¬ 
tion for the dipole. Conversely, if unquestionable cosmological observations 
conclude that our local motion is in disagreement with the CMBR dipole, 
then it would suggest that we live in a flat universe and that ultra large scale 
isocurvature perturbations must exist. 
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